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S : ABSTRACT 

^ ' We present the evolution of the radial expansion of SNR 1987A as measured 

using Chandra X-ray observations taken over the last 10 years. To characterize 
the complex structure of the remnant and isolate the expansion measurement, 
we fit the images to several empirical models including: a simple circular torus, 
W ' a torus with bilateral lobes, and a torus with four tangentially extended lobes. 

K ■ We discuss the results of this measure in the context of the overall evolution of 

• the supernova remnant, for which we believe we have measured the end of the 

^ ' free expansion phase and its transition to the adiabatic phase (at least along the 

i-{ • equatorial ring). The timing of this event is in agreement with early predictions 

52 ■ of the remnant evolution. 

> 

■ 1. Introduction 

0^ ■ 

^ . Supernova 1987A is the nearest and consequently the most extensively studied supernova 

QQ . (SN) in modern astronomy. While it has led to important insights into progenitor systems 

^ . and supernova evolution, it is by no means a typical ch aracteristic event. Th e progenitor 

O ■ was seen in archival plates as blue supergiant Sk -69 202 ( 1 White fc Malinlll987l ) rather than 
the expected red supergiant progenitor of a Type II SN. The extent of the stellar winds from 

^ , both the red and blue supergiant phases suggest that the sys tem transitioned from a red 



supergiant sometime in the last ~ 40,000 years (lArnett et al.l 119891 ). This transition may 



explain some of the unusual properties of the supernova remnant (SNR). 

The ring-like structure visible and evolving over time in X-ray, optical, and radio ob- 
servations was first discovered to have "clumpy" structure using the high-resolution Hubble 
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Space Telescope images (jPlait et al.lll995l ). The many spots seen in the optical appear to cor- 
respond to the less resolved soft X-ray spots. The X-ray and radio image struc tures are much 



simp ler than the complex triple-ring structure seen in optical observations (IBurrows et al. 



19951 ): the X-rays and radio emission are only detected from the central inner ring. This 
ring appears to be a structure of hot gas that is flattened into a round disk. It appears 
elliptical due to its inclination, is gradually brightening in all bands, shows increasing struc- 
tural details, and is expanding radially. The optical emitting region of the inner ring is 
slightly outside that of the X-ray and radio emitting regions. The soft X-ray emission from 
SNR 1987A is thought to be produced by shocked hot gas that lies between the forward 
and reverse shocks of the SN blast wave, with the localized spots of emission around the 
inner ring being du e to regions of enhanced density protrud ing from the inner surface of the 
circumstellar ring (jPark et al.l 120021 : ISugerman et al.l |2002| ) . The inner and outer rings are 
thought to be produced by asy mmetric stellar winds shed at different stages in the evolution 
of the massive progenitor star (ILuo fc McCraylll99ll : IChevalier fc Dwarkadaslll995l ). 



The soft X-ray emission has been increasing in intensity as a result of the blast wave 
shocking and continually sweeping up the surrounding circumstellar material. The structure 
of the remnant has evolved from emission largely contained in these small regions to overall 
brightening of the ring during the progressive observations of the resolved remnant. As the 
blast wave lights up the dense inner r ing, the X-ray spe ctrum may become dominated by 
the density structure of the ring itself (jPark et al.ll2006al ). Figure [1] shows the 0.3 — 8.0 keV 
Chandra images at each epoch (Tabled]), showing how the X-ray emission has transitioned 
in recent years from being dominated by the shocked gas in t he dense hot s pots to being 
dominated by the gas between the forward and reverse shocks (jMcCrayll2003l ). 



The ROSAT X-ray light curve increased linearly for approximately 10 years after the 
SN explosion. Later Chandra observations indicated that there was a deviation beginning 
in ~ 1997 (~ 3700 days after the SN) from the ROSAT linear extrapolation of the flux, 
which results in a nearly exponential increase in the soft X-ray light curve (Figu re [21). Th is 



transition was temporally co incident with the appearance of the first optical spot (jPark et al. 



2OO2I : ISonneborn et al.lll998l ) . This rapid flux increase was followed by an upturn in the X-ray 
light curve in 2004 (~ 6200 days after the SN), suggesting that the blast wave had reached 
the main body of the inner ring. This brightening was accompanied by substantial softening 
of the overall X-ray spectrum, brightening of hot spots in the ring, filling in of the X-ray 
ring, and a decrease in radial expansion velocity of the ring. 

In this paper we describe a new method for measuring the size of the remnant from the 
X-ray images, and use that to put it in the context of the structure and overall evolution 
of the system. We describe in ^ the predictions in the literature prior to our Chandra 
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campaign, the observations in ^ our radial expansion measurement and image models in 
§11 the results in ^ the discussion in §6l and our conclusions in ^ 



Predictions 



Early observations of the circumstellar ring and measurements of the ejecta velocities 
prorapted predictions for w hen the blast wave would encounter the inner ring (jChevalier &: Liang 
19891 : ILuo fc McCraylll99ll ). This encounter was predicted to occur between 8 and 16 years 
after the SN and would result in dramatic re-brightening as a result of the strong shocks 
produced by the blast wave interact i ng wi th the inner ring. This prediction was later 
revised when I Chevalier fc Dwarkadad (119951 ) suggested that an H II region may exist in- 
side the inner circumstellar ring. It is needed to reconcile the expected low mass loss rate 
(M ~ 7.5 X lO~^M0yr~^) of the blue supergiant progenitor with the early radio free-free 
absorption and the X-ray emission levels. The H II region was probably formed in material 
swept up by the red supergiant wind that was ionized by the blue supergiant star. This has 
the effect of delaying the date of the encounter with the inner ring to year 2005 ± 3 because 
the blast wave is slowed as it swept up this additional material. 



Borkowski et al.l ( 119971 ) modeled the SN blast wave impact with the inner circumstellar 



ring and the resulting reverse shock. They predicted that the impact would be characterized 
by a gradual increase in X-ray luminosity. As the shock front engulfed the inner ring, it was 
expected to gradually brighten by as much as 3 orders of magnitude in the optical, UV, and 
X-ray. The X-ray spectrum would also change depending on the detailed shock structure of 
the interaction with the inner ring. The overall soft X-ray flux evolution is consistent with 
these predictions, showing a steep increase in the light curve since ~ 2000 (see Figure [2]). 
However, these predictions were made before the emergence of the first optical spot, and 
therefore could not have accounted for how the higher density protrusions would affect the 
observational behavior. 

After ~ 1997 (~ 3700 days after the SN) the soft X-ray light curve began to deviate from 
the linear increase observed for the prior 10 years. This brightening was initially confined to 
the localized hot spots evident in HST images in the northeastern and southeastern regions 
of the inner ring. As more optical hot spots appeared, the X-ray emission also spread 
around the ring and the X-ray luminosity increased dramatically. The brightening of the 
entire ring of SNR 1987A appears to have begun in early 200 4 and the events predicted 
are underway. This is evident in the brightenin g in X-rays (jPark et al.l |2005| ) , infrared 
( iBouchet et al.l 120061). and optical (|McCravl 120051) . the continuous softening of the X-ray 
spectrum (jPark et al.ll2006bl : IZhekov et al.ll2005l ). deceleration of the radial expansion of the 
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SNR (Park et al 



( iMichael et al 



2004 



2002 



Dewey et al 



2008 



Zhekov et al.ll2005r). and slower shoc k velocities from line broadening 



Zhekov et all 120091 ) . 



3. Chandra Observations and Data Reduction 

We have observed SNR 1987A with Chandra on 23 occasions between 1999 October 
and 2009 January, and our monitoring observations continue at approximately 6 month in- 
tervals. Eighteen of the observations were carried out using the Advanced CCD Imaging 
Spectrometer (ACIS) back-illuminated S3 CCD. ACIS alone provides high-resolution images 
with moderate non-dispersive energy resolution. When used in conjunction with one of the 
gratings, either the High Energy Transmission Grating (HETG) or the Low Energy Transmis- 
sion Grating (LETG), ACIS provides high quality spectra and an undispersed zeroth-order 
image with moderate energy resolution. Six of our obse rvations used the gratings: HETG 



1999 October and 2007 April, 2008 July, 2009 January, (jBurrows et al.ll2000l : Michael et al. 



20021: 



20051. 



Dewev et al.ll2008h ) and LETG (2004 September and 2007 September, ( 



200f 



Zhekov et al. 



^, 2OO9I )). In the standard LETG configuration, the long- wavelength oxygen lines 



are positioned on the sensitive ACIS-S3 chip, which consequently places the zeroth order 
image 1.5' off-axis, blurring it by Chandra^s off-axis Point Spread Function (PSF). This 
blurring degrades the image, and we therefore exclude the LETG observations from our 
measurement of the radial expansion rate. The HETG, however, does not suffer from this 
effect. The differences in the images between the bare- ACIS and the zeroth-order HETG 
have not been calibrated, though no significant systematic effects (as seen with the LETG) 
are expected. The first HETG observation is also our first monitoring observation, therefore 
it is an important reference point for our analysis. The deep HETG observation taken in 
2007 was not part of our regular monitoring observations and was performed in between 
regularly scheduled monitoring observations. As the remnant brightens, we have recently 
taken calibration observations comparing the bare- ACIS and HETG zeroth-order images (in 
2008-7 and 2009-1). In order to avoid photon pipe-up, we will transition future observations 
to HETG only. These closely spaced HETG observations to monitoring ACIS observations 
revealed consistent radii between the ACIS-only and HETG measurements. They do not 
add new information to this study and therefore do not have a strong affect on the radial 
expansion velocity fits. We exclude all gratings observations except the 1999 October HETG 
observation from the rest of radial expansion analysis presented in this paper. The details 
of the 19 Chandra observations used in this work are listed in Table [TJ 



The technique u s ed to create our images is the same as was used in our previous works 
JSurrows et allboool : [Park et al.lbood l2004j . boOGal lbll The event data used to create the 
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images were filtered to use the energy range 0.3 — 8 keV, and ASCA grades 02346. Flaring 
pixels were removed. We corrected for the charge transfer inefficiency (CTI) using the 
methods developed by iTownsley et al.l (120001 ). The small angular extent of the remnant 
(~ 1'.'6) makes any image analysis difficult even with the superb spatial resolution of Chandra 
(ACIS CCD pixel size ~ 0M9 2). To better s ampl e the image, we appli ed the sub-pixel 
resolution method developed by lTsunemi et al.l (120011 ) and lMori et al.l (1200 ll ). This technique 
uses the feature that information can be obtained on the location within the pixel where the 
photon was absorbed by examining the distribution of the events split between pixels. The 
morphology of the splitting can be traced to a centroid within the pixel for events spread 
across more than 2 pixels. This technique improves the angular resolution of the images by 
~ 10%. 

Several other factors also help to improve the image resolution with the use of the sub- 
pixel resolution technique. The choice of using one of the ACIS back-illuminated CCDs 
leads to a larger fraction of split-events. The back illuminated chips are also less affected by 
the poor charge transfer inefficiency fo r observations obtained after the particle damage in 
September 1999 (iTownsley et al.ll2000l ). Another source of better resolution comes from the 
intentional dithering motion of the observatory which causes a target image to move across 
the CCD surface forming a Lissajous pattern. This leads to more opportunities of forming 
split-pixel events in different areas of the image. Standard processing adds a randomization 
of event positions within each pixel to prevent marked pixelization of the images. We turn 
off this randomization in the data processing to prevent any loss of positional information. 

To further improve the effective a ngular resolution of the i mage s, we deconvolve them 
with a maximum likelihood algorithm (IRichardsonl 1 19721 : iLucyl Il974l ) using the on-axis de- 
tector PSF. We used a scale of 0'.'125 sky pixels and then smoothed the images by convolving 
with a Gaussian (~ O'.'l FWHM). The resulting 0.3 — 8.0 keV broadband images of SNR 
1987A from each of the epochs of our observations are shown in Figure [U 



4. Radial Expansion 



Previo us works to model the X-ray radial expansion rate of SNR 1987A (jPark et al. 



2OO2I . I2OOJ) used a simple Gaussian model to fit the radially averaged profile. This approach 



yielded an expansion velocity estimate of 4000 — 5000 km s~^. The utility of this simple 
method was limited by the loss of information on the detailed features such as the lobes, 
the overall asymmetric intensity, as well as ellipticity and inclination of the inner ring. To 
more accurately quantify the radial expansion, we developed a more realistic image model 
to estimate the SNR's radius. 
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The goal for our empirical models is to characterize the Chandra X-r ay images in a 



more realistic way and in a similar spirit as tho se of the radio studies (iGaensler et al. 



19971 : [Manchester et al.l |2002| . l2005l : iNg et al.l |2008| ) in order to obtain a measurement for 



the radial expansion of th e remnant. The inner ring is inclined by 43° to the line of sight 
( ICrotts fc Heathcotelll99ll ). The bulk of the X-ray emission originates from the "disk" con- 
taining the inner ring, rather than from a spherical volume (jZhekov et al.l 120051 ) . The X-ray 
morphology of SNR 1987A appears to be an elliptical ring with 3-4 enhanced lobes (Figure 
[T]). The real number of X-ray spots is likely larger, but the Chandra spatial resolution pro- 
hibits the distinguishing of finer image details. To remove the apparent ellipticity prior to 
fitting the observed images with our models, we deproject the images by the 43° inclination. 
After the deprojection, the images are approximately circular (e.g.. Figure [3^ andlHb). We 
fit three simple empirical models to the image data in an attempt to measure the radial ex- 
pansion velocity, with the latter two models removing the biases created by the brightening 
hot spots. 



4.1. Torus Only Model 



The torus model is the simplest of the three models. It assumes that the images can be 
approximated by a circularly symmetric smooth ring with a peak radius of ro and a Gaussian 
width of CTj-Q, with the following form: 



r = [{x- XoY + iy- VoY] 



Mi 



torus 



No + Nr exp 



1/2 



2al 



(2) 



where x and y are the coordinates of each pixel in the model image, xq and yo are the 
offsets of the ring center from the image center, Nq is the background image level, A^^ is the 
normalization, and Mt^rus is the resulting model image. An example of a fit to this model is 
shown in Figure [St. 



4.2. Bilateral Lobes -|- Torus Model 



This model consists of the torus as described in §4. II plus 2 bilatera l lobes. The addition 
of bilateral lobes approximates the morpholo gy of the radio remnant (IGaensler et al.l 119971 : 
Manchester et al.l 120021 . l2005l : iNg et al.ll2008l ). The morphology of the X-ray remnant in the 
earliest epochs resembled that of the radio remnant with emission approximated by opposite 
bilateral regions of increased brightness. As the X-ray spots have brightened, additional 
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lobes are required to adequately model the images. Therefore, this model works best in 
the early epochs and poorly in late ones. We required the lobes to be positioned with one 
on each the East and West side of the remnant and allow them to vary in intensity (Nl-), 
position angle (6'j), and azimuthal width (o"t-). The torus is also allowed to vary in xq, z/q, tq 
and (Jro- The lobes are then multiplied by the torus which puts them at a radius of ro and 
a radial width of aro, taking the form: 



L, exp 



where i represents each individual lobe. We combine the torus and the lobes as: 

M2lobes = Mtorusi^ + Li + L2) 

An example of the fit to this model is shown in Figure [3ll. 



(3) 



(4) 



4.3. Quadrilateral Lobes + Torus Model 

A comparison of Figure [3)d and [3li (bilateral model) clearly shows that the bilateral 
model fails to reproduce the 4 major hot spots seen in the X-ray image. Because there 
appears to be one such hot spot in each quadrant of the image, we next attempt to model 
the data using 4 tangentially extended lobes superposed on the torus described in 14.11 with 
one lobe in each image quadrant (NW, NE, SE, SW). However they are allowed to vary in 
intensity (Nt.), position angle {9i) within the quadrant, and azimuthal width (0"^.), where i 
represents the properties of each individual lobes. The torus is also allowed to vary in xq, 
Hq, tq and cr^o- The lobes are then multiplied by the torus which puts them at a radius of 
ro and a radial width of aro, taking the forms of Equation [31 We combine the torus and the 
lobes as: 

M^lobes = Mtorusil + Li + L2 + L3 + L4) (5) 

An example of the fit to this model is shown in Figure [Sf. 



5. Results 

We fit the deprojected images from each epoch to each of the three models using a 
Levenberg-Marquardt two-dimensional least-squares IDL fit procedure. Our models are able 
to sufficiently represent the data in all cases, with improved reduced model complexity 
increases. Example image fits to all models for sample observation 11 are presented in Figure 

El 
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The main goal of this study is to measure the expansion rate of the remnant and follow its 
evolution. We use the parameter tq measured in each model fit to characterize this behavior. 
We convert rn into physical velocity units using the distance to SNR 1987A (51.4 ± 1.2 kpc, 
Panagial l2003l ). The expansion measures do not increase linearly throughout, but rather 
are better fit by a linear fit that extends until ~ 6000 days after the SN, and then turns 
over to a shallower slope (Table [21 Figure H]). Table [3] lists the resulting radial expansion 
fit parameters for each set of the image model fits. All three models provide consistent fits 
to the radial expansion measurements with the tightest constrains on tq from the model 
that best represents the images (quadrilateral lobes plus torus). The results of these fits are 
shown Figure m 

Our estimates of the radial exp ansion rates before day 6000 are ~ 65% larger than 
previous estimates (jPark et al.ll2004l ). This is most likely due to the two different methods 
used to measure the radius of the SNR. Our radius estimates use more sophisticated image 
model fits than the simple Gaussian radial profile of the X-ray ring used by the previous 
work, and thus should be more realistic and reliable. We also considered the 43° inclination 
effect of the inner ring causing it to appear elliptical, whereas our previous work assumed 
a simple circular geometry. For comparison, if we re-incline the measured radial expansion 
rate, the corresponding expansion velocity is ~ 5500l}^oQkms~^ prior to day ~ 6000 and ~ 
1200: 



-460 
-440 



kms~^ afterwards. This is in plausible agreement with the previous measurements. 
Hydrodynamic models indicate a forward shock velocity of ~ 3500 — 4500 km s~^ until day ~ 
6000, and then a dramatic deceleration down to ~ 1000 — 1500 km s~^ afterwards (iDwarkadas 



20071 ). This overall change of the shock velocity is also in good agreement with the results 
from our expansion rate estimates. 



The most recent radio radial expansion measurements by iNg et al.l (120081 ) show a con- 
stant remnant expansion rate of 4000 ± 400 km s^^ over a range of 1992-2008. However, the 
radio remnant radius is systematically larger than the X-ray remnant radius by ~ 10 — 20%. 
Interestingly, after day ~ 5400, the radio rern nant radius is larger than the radius of the 
optical inner ring (0'.'83, ISugerman et al.ll2002l ). This discrepancy in the SNR size between 
the X-ray an d radio data is largely caused by the different image modeling techniques as 
discussed by iNg et al.l (120081 ) . It may also indicate that the X-ray and radio emission are 
produced in different physical regions. The radio emission is likely due to synchrotron pro- 
cesses in between the forward and reverse shocks (IManchester et al.ll2005l ). while the X-rays 
are dominated by thermal emiss ion from the trans mitted and reflected shocks produced by 
the shock-inner ring interaction ( Zhekov et al.|[2009l ). The 3- dimensio nal effects and emission 
regions also likely differ between these components. iNg et al" (2008) predict that they will 
soon see a deceleration in the radio remnant expansion similar to that we have observed 
from the X-ray remnant due to the forward shock interacting with the optical inner ring. 
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6. Exploring Methods with Simulated Images 

In order to test the reliability of our image modeling and radial expansion measures, we 
created a series of simulated images using the Model of AXAF Response to X-rays (MARX) 
tool. The simulated images were meant to approximate the real images obtained at days 
~ 4600, 5200, 6200, and 7200 days after the SN explosion, assuming the estimated radius 
for each epoch. They were generated with different observed inclinations and properties, 
specifically degrees inclination (face-on) , and 45 degrees inclination both with and without 
4 spots. We applied the Chandra PSF to the simulated images to mimic observing effects, 
and ran them through the same image deconvolution process as done with the real data. 

We processed these simulated images through the same algorithms used to fit the ob- 
served images for all of our image models. The results showed a systematic offset in the 
best-fit radius, which is ~ 10% smaller than the assumed input values, yet gave an accurate 
measure of the radial expansion rate within measurement errors. We attribute this system- 
atic offset to differences in the treatment of the torus width (c-r,,). The simulated images were 
constructed assuming that the torus and any changes to the width over time are not spatially 
resolved. Applying this same assumption to the simulation fits, we froze the torus width to 
approximately the value of the effective PSF for the deconvolved images (~ 0'.'2 — 0'.'3), and 
redid the fits. This effectively ehminated the systematic offset and reproduced the input 
radii accurately. 

In the application to the real data, when leaving the torus width as a free parameter, 
the model fits may underestimate the radius of the real remnant images by as much as 
10%. However, the free torus width fits appear to reveal a widening torus with values at a 
similar size to the PSF. Any real physical changes to the torus width would be ignored if the 
parameter is frozen to a set value. The width is not independent of the radii measurement 
that we are inherently trying to measure. The fits to the real data are statistically improved 
if we assume that the torus has some real width that may change over time and the width 
parameter is left free to vary. The fitted torus width (Gaussian half-width) ranges from 
~ 0'.'12 — 0'.'21 which is significantly smaller than the ring radius and comparable to the 
effective PSF of the deconvolved Chandra images. 

We also assume a 2-dimensional structure in all of our model fitting even with the 
deprojecting correction. In reality, X-ray emission most likely originates from a 3-dimensional 
structure with the shock component being perpendicular to the disk containing the inner 
ring, making contributions to the torus width unknown. The real 3-dimensional geometry 
of the unresolved inner ring is not yet well understood. The potential ~ 10% radii offset is 
systematic and may consequently have the same effect on the expansion velocities. However, 
the measurement errors are larger than this systematic offset. 
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The first two images are most strongly affected by freezing the torus width because 
they suffer from low count statistics (lower by an order of magnitude than those at other 
epochs), and contain the lowest contributions of the X-ray emission from the shock-inner 
ring interaction. Therefore, this 3-dimensional effect could have been relatively large in 
these early images. To test this effect, we fit the radial expansion rate using the frozen 
torus width with and without the first two images, with the latter analysis leading to a more 
accurate reproduction of the results from the free torus width fits. When including the first 
two observations with a frozen torus radius {tq ~ 0.15"), the initial expansion velocity is 
~ 12000 km s~^, fiattening to a velocity of ~ 1400 km s~^ at ~ 5800 days. The expansion 
measurement excluding the first two observations yields an initial expansion velocity of 
~ 6000 km s~^ fiattening to a velocity of ~ 1700 km s~^ at ~ 6200 days. Therefore, it is 
these early images which are infiuencing the fits and are least conducive to a frozen torus 
width. Even with a free torus width, the first two images are the noisiest and provide the 
poorest fits. The first two images have some effect on the accuracy of the expansion rate, 
but do not change the overall conclusion that we observe a slowing of the expansion rate 
starting at ~ 6000 days. 



7. Conclusions 



Mult i- wavelength studies at frequent epochs are carefully capturing the evolution and 
progression of SNR 1987A. As the evolution continues, it reveals new regions of emission 
and probes a larger volume of the circumstellar environment. We believe that the pre- 
dicted events (IChevalier fc Liangj|l989l : IChevalier fc Dwarkadaslll995l : ILuo fc McCraylll991 



Borkowski et al.lll997l ) of the SN blast wave engulfing the inner circumstellar ring are occur- 
ring. The earlier optical and soft X-ray emission spots were a result of higher density gas 
protruding inward from the main ring. In recent years, we have seen the entire ring begin 
to light up. The current optical and X-ray emission is a composite of the shocks heating the 
higher and lower density regions of the clumpy ring structure. 

Our measurement of the deceleration of the radial expansion of SNR 1987A from 
~ 8000 km s~^ to ~ 1600 km s~^ at ~ 6000 days along with the upturn of the X-ray lu- 
minosity, fractional contribution from the decelerated shock to the soft X-ray fiux, and the 
low shock velocities from line broadening are all pointing to indications of the SN blast wave 
encountering the main body of the circumstellar ring around day 6200. The evidence of the 
blast wave encountering the inner ring are already apparent in the X-ray observations, and 
we are yet to see the strong indications from the radio or optical observations. The details 
of the different emission regions and their interactions with the different shock components 
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will become clearer as continued optical, X-ray, and radio observations probe the remnant. 

Continued monitoring of SNR 1987A with Chandra will reveal more details in the struc- 
ture of the circumstellar environment around the remnant. We expect to see additional 
brightening as the blast wave shock front heats the material in the ring and as the reverse 
shock region continues to grow in volume. As the remnant expands and brightens, we hope 
to further resolve structures and to understand the inhomogeneous nature of the circumstel- 
lar environment. We expect the expansion to continue to slow as it sweeps up more material 
and transitions from the free expansion phase to the adiabatic expansion phase once it has 
swept up a mass of material comparable to the SN ejecta. Uncertainties in this ejecta mass 
make predictions for this time frame difficult. However, SNR 1987A has an unusually high 
density environment as evident by the effects on the original SN explosion. Therefore this 
transition may occur earlier than the typical ~ 200 year timescale. 

We may also expect to see additional deceleration of the radial expansion while the blast 
wave sweeps up more of the circumstellar material. Only the material that was photo-ionized 
by the UV flash from the SN is visible in the optical i nner ring, and there is evidence suggest- 



ing th at there may be more mass that was not ionized (jCrotts et al.lll989l : lEnsman fc Burrows 



I992I ). The inner ring has an ionized mass of ~ O.O4M0, which is probably only the inner 



surface of a larger gas complex. As the hot spots brighten and produce ionizing r adiation. 



event ually the undetected outer material will become visible as an emission nebula (iMcCray 



20031 ). The UV and optical light echoes seen several months after the SN event hinted at 
the structure of the dust in the outer circumstellar environment. It will be interesting to 
continue to use the radial expansion measurement to probe the density structure of the ring. 



This work is supported by Chandra SAO grants GO7-8069X and GO8-9076X. S.P. was 
also supported in part by SAO grant SV4-74018. 



- 12 - 



REFERENCES 

Arnett, W. D., Bahcall, J. N., Kirshner, R. P., & Woosley, S. E. 1989, ARA&A, 27, 629 
Borkowski, K. J., Blondin, J. M., & McCray, R. 1997, ApJ, 477, 281 

Bouchet, R, Dwek, E., Danziger, J., Arendt, R. G., De Buizer, I. J. M., Park, S., Suntzeff, 
N. B., Kirshner, R. P., & Challis, P. 2006, ApJ, 650, 212 

Burrows, C. J., Krist, J., Hester, J. J., Sahai, R., Trauger, J. T., Stapelfeldt, K. R., Gallagher, 
III, J. S., Ballester, G. E., Casertano, S., Clarke, J. T., Crisp, D., Evans, R. W., 
Griffiths, R. E., Hoessel, J. G., Holtzman, J. A., Mould, J. R., Scowen, P. A., Watson, 
A. M., & Westphal, J. A. 1995, ApJ, 452, 680 

Burrows, D. N., Michael, E., Hwang, U., McCray, R., Chevalier, R. A., Petre, R., Garmire, 
G. P., Hoh, S. S., & Nousek, J. A. 2000, ApJ, 543, L149 

Chevaher, R. A., & Dwarkadas, V. V. 1995, ApJ, 452, L45 

Chevalier, R. A., & Liang, E. P. 1989, ApJ, 344, 332 

Crotts, A. P., & Heathcote, S. R. 1991, Nature, 350, 683 

Crotts, A. P. S., Kunkel, W. E., & McCarthy, P. J. 1989, ApJ, 347, L61 

Dewey, D., Zhekov, S. A., McCray, R., & Canizares, C. R. 2008, ApJ, 676, L131 

Dwarkadas, V. V. 2007, in American Institute of Physics Conference Series, Vol. 937, Amer- 
ican Institute of Physics Conference Series, ed. S. Immler & R. McCray, 120-124 

Ensman, L., & Burrows, A. 1992, ApJ, 393, 742 

Gaensler, B. M., Manchester, R. N., Staveley-Smith, L., Tzioumis, A. K., Reynolds, J. E., 
& Kesteven, M. J. 1997, ApJ, 479, 845 

Lucy, L. B. 1974, AJ, 79, 745 

Luo, D., & McCray, R. 1991, ApJ, 372, 194 

Manchester, R. N., Gaensler, B. M., Staveley-Smith, L., Kesteven, M. J., & Tzioumis, A. K. 
2005, ApJ, 628, L131 

Manchester, R. N., Gaensler, B. M., Wheaton, V. C, Staveley-Smith, L., Tzioumis, A. K., 
Bizunok, N. S., Kesteven, M. J., & Reynolds, J. E. 2002, Publications of the Astro- 
nomical Society of Austraha, 19, 207 



-13- 



McCray, R. 2003, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 598, Supernovae 
and Gamma-Ray Bursters, ed. K. Weiler, 219-240 

McCray, R. 2005, in lAU CoUoq. 192: Cosmic Explosions, On the 10th Anniversary of 
SN1993J, ed. J.-M. Marcaide & K. W. Weiler, 77 

Michael, E., Zhekov, S., McCray, R., Hwang, U., Burrows, D. N., Park, S., Garmire, G. P., 
Holt, S. S., & Hasinger, G. 2002, ApJ, 574, 166 

Mori, K., Tsunemi, H., Miyata, E., Baluta, C. J., Burrows, D. N., Garmire, G. P., & 
Chartas, G. 2001, in Astronomical Society of the Pacific Conference Series, Vol. 251, 
New Century of X-ray Astronomy, ed. H. Inoue & H. Kunieda, 576 

Ng, C.-Y., Gaensler, B. M., Staveley-Smith, L., Manchester, R. N., Kesteven, M. J., Ball, 
L., & Tzioumis, A. K. 2008, ApJ, 684, 481 

Panagia, N. 2003, ArXiv Astrophysics e-prints 

Park, S., Burrows, D. N., Garmire, G. P., Nousek, J. A., McCray, R., Michael, E., & Zhekov, 
S. 2002, ApJ, 567, 314 

Park, S., Zhekov, S. A., Burrows, D. N., Garmire, G. P., & McCray, R. 2004, ApJ, 610, 275 

Park, S., Zhekov, S. A., Burrows, D. N., Garmire, G. P., Racusin, J. L., & McCray, R. 2006a, 
ApJ, 646, 1001 

Park, S., Zhekov, S. A., Burrows, D. N., & McCray, R. 2005, ApJ, 634, L73 

Park, S., Zhekov, S. A., Burrows, D. N., Racusin, J. L., McCray, R., & Borkowski, K. J. 
2006b, in ESA Special Publication, Vol. 604, The X-ray Universe 2005, ed. A. Wilson, 
335-340 

Plait, P. C, Lundqvist, P., Chevaher, R. A., & Kirshner, R. P. 1995, ApJ, 439, 730 

Richardson, W. H. 1972, Journal of the Optical Society of America (1917-1983), 62, 55 

Sonneborn, G., Pun, C. S. J., Kimble, R. A., Gull, T. R., Lundqvist, P., McCray, R., Plait, 
P., Boggess, A., Bowers, C. W., Danks, A. C, Grady, J., Heap, S. R., Kraemer, S., 
Lindler, D., Loiacono, J., Maran, S. P., Moos, H. W., & Woodgate, B. E. 1998, ApJ, 
492, L139 

Sugerman, B. E. K., Lawrence, S. S., Crotts, A. P. S., Bouchet, P., & Heathcote, S. R. 2002, 
ApJ, 572, 209 



-14- 



Townsley, L. K., Broos, P. S., Garmire, G. P., & Nousek, J. A. 2000, ApJ, 534, L139 

Tsunemi, H., Mori, K., Miyata, E., Baluta, C, Burrows, D. N., Garmire, G. P., & Chartas, 
G. 2001, ApJ, 554, 496 

White, G. L., & Malin, D. F. 1987, Nature, 327, 36 

Zhekov, S. A., McCray, R., Borkowski, K. J., Burrows, D. N., & Park, S. 2005, ApJ, 628, 
L127 

— . 2006, ApJ, 645, 293 

Zhekov, S. A., McCray, R., Dewey, D., Canizares, C. R., Borkowski, K. J., Burrows, D. N., 
& Park, S. 2009, ApJ, 692, 1190 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 15 - 



Table 1. Chandra Observations of SNR 1987A 



Observation # 


Obs Date 


Age* 


Instrument* 


Exposure 
(ks) 


Source 
Counts* 


1 


1999 October 6 


4608 


ACIS-S + HETG^ (3.1 s) 


116.1 


690 


2 


2000 January 17 


4711 


ACIb-b3 (3.2 s) 


8.6 


607 


3 


2000 December 7 


5036 


ACIb-b3 (3.2 s) 


98.8 


9030 


4 


2001 April 25 


5175 


ACIb-b3 (3.2 s) 


17.8 


1 cir\r\ 

1800 


5 


2001 December 12 


5406 


ACIS-S3 (3.1 s) 


49.4 


6226 


6 


2002 May 5 


5560 


AC1S-S3 (3.1 s) 


44.3 


6427 


7 


r\r\r\r\ T~\ i o "i 

2002 December 31 


5790 


AC1S-S3 (3.1 s) 


49.0 


9277 


8 


2003 July 8 


5979 


AC1S-S3 (3.1 s) 


45.3 


9668 


9 


2004 January 2 


6157 


ACIS-S3 (3.1 s) 


46.5 


11856 


10 


2004 July 22 


6359 


ACIS-S3 (1.5 s) 


48.8 


17979 


11 


2005 January 12 


6532 


ACIS-S3 (0.4 s) 


48.3 


24939 


12 


2005 July 14 


6716 


ACIS-S3 (0.4 s) 


44.1 


27048 


13 


2006 January 28 


6914 


AC1S-S3 (0.4 s) 


42.3 


30940 


14 


2006 July 28 


7094 


AC1S-S3 (0.4 s) 


36.4 


30870 


15 


2007 January 19 


7270 


AC1S-S3 (0.4 s) 


33.5 


32798 


16 


2007 July 13 


7445 


ACIS-S3 (0.4 s) 


25.7 


27945 


17 


2008 January 10 


7626 


ACIS-S3 (0.2 s) 


9.3 


12008 


18 


2008 July 4 


7802 


ACIS-S3 (0.2 s) 


8.6 


12119 


19 


2009 January 5 


7987 


ACIS-S3 (0.2 s) 


6.0 


9204 



*Days since SN explosion 

*Zeroth order data from grating observations 

*0.3 - 8.0 keV band 

*Time in parentheses is the ACTS frame-time. 
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Tablc 2. Model Parameter Fits 



Age* 


Torus Only 


Torus±2 lobes 


Torus±4 lobes 








rn V ^ij^^ / 






/ (J \ (All. 1 




0.12 ±0.01 


0.60 ±0.02 


0.12 ±0.01 


0.60 ±0.02 


0.12 ±0.01 


0.60 ± 0.02 


471 1 


0.14 ±0.01 


0.59 ±0.02 


0.14 ±0.01 


0.59 ±0.02 


0.14 ±0.01 


0.59 ±0.02 


(JUOU 


0.19 ±0.01 


0.65 ±0.01 


0.18 ±0.01 


0.66 ±0.01 


0.18 ±0.01 


0.66 ±0.01 


■^1 7^ 

OJ- i o 


0.17 ±0.01 


0.66 ± 0.02 


0.16 ±0.01 


0.64 ± 0.02 


0.17 ±0.01 


0.67 ± 0.02 


5406 


0.19 ±0.01 


0.68 ± 0.01 


0.19 ±0.01 


0.68 ± 0.01 


0.19 ±0.01 


0.69 ± 0.01 


5560 


0.19 ±0.01 


0.68 ±0.01 


0.19 ±0.01 


0.67 ±0.01 


0.19 ±0.01 


0.69 ±0.01 


5790 


0.22 ±0.01 


0.69 ±0.01 


0.22 ±0.01 


0.71 ±0.01 


0.22 ±0.01 


0.72 ±0.01 


5979 


0.20 ± 0.01 


0.73 ±0.01 


0.21 ±0.01 


0.72 ± 0.01 


0.21 ±0.01 


0.73 ±0.01 


6157 


0.22 ±0.01 


0.73 ± 0.01 


0.22 ±0.01 


0.73 ±0.01 


0.22 ±0.01 


0.74 ± 0.01 


6359 


0.22 ±0.01 


0.74 ± 0.01 


0.22 ± 0.01 


0.74 ± 0.01 


0.22 ± 0.01 


0.74 ± 0.01 


6532 


0.22 ±0.01 


0.74 ± 0.01 


0.21 ±0.01 


0.74 ± 0.01 


0.21 ±0.01 


0.74 ± 0.01 


6716 


0.21 ±0.01 


0.75 ±0.01 


0.21 ±0.01 


0.74 ±0.01 


0.21 ±0.01 


0.75 ±0.01 


6914 


0.22 ±0.01 


0.74 ±0.01 


0.22 ±0.01 


0.74 ±0.01 


0.22 ±0.01 


0.75 ±0.01 


7094 


0.21 ±0.01 


0.76 ±0.01 


0.21 ±0.01 


0.75 ±0.01 


0.21 ±0.01 


0.76 ±0.01 


7270 


0.21 ±0.01 


0.76 ± 0.01 


0.21 ±0.01 


0.76 ± 0.01 


0.21 ±0.01 


0.76 ± 0.01 


7445 


0.20 ±0.01 


0.77 ± 0.01 


0.21 ±0.01 


0.77 ± 0.01 


0.21 ±0.01 


0.77 ±0.01 


7626 


0.21 ±0.01 


0.76 ± 0.01 


0.21 ±0.01 


0.76 ± 0.01 


0.21 ±0.01 


0.76 ± 0.01 


7802 


0.20 ±0.01 


0.77 ±0.01 


0.20 ± 0.01 


0.76 ±0.01 


0.20 ±0.01 


0.76 ± 0.01 


7987 


0.20 ±0.01 


0.78 ± 0.01 


0.20 ± 0.01 


0.78 ±0.01 


0.20 ±0.01 


0.78 ± 0.01 



*Days since SN explosion 

Note. — Errors are la confidence intervals. 



-17- 



Table 3. Radial Expansion Velocity Fits 





Va 

(km s"-*-) 


turn over time 
(days) 


Vb 

(km s"-*-) 




Torus only 
Torus + 2 Lobes 
Torus + 4 Lobes 


6908+}^83 
7Q47+1460 

' "J^' -1204 
7COQ+2139 


6135_284 


lOfi 1+602 

2037+^^^ 


0.70 
0.53 
0.54 



Note. — Va and vt, are calculated using the slopes of the broken 
linear fit with 15 DOF. All errors are 90% confidence. 
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Fig. 1. — Broadband (0.3 — 8.0 keV) images of SNR 1987A from each Chandra observation 
taken on date (YYYY-MM) as labeled. The images are all the same spatial scale with slight 
uncertainty in rotations and centers due to unknown absolute astrometry on these scales. 
The images are scaled arbitrarily to show the morphological details of the evolving remnant. 
The flux has actually increased dramatically over the course of this monitoring campaign, 
as shown in Figure [21 
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Fig. 2. — Soft X-ray (0.5 — 2 kev) light curve from ROSAT and Chandra observations. 
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Fig. 3. 
image, 
model, 



, — Images and fits to models for observation 11. a: original deconvolved sub-pixel 
b: deprojected image, c: fit to torus only model, d: fit to the torus plus two-lobes 
e: fit to torus plus four-lobes model. 
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Fig. 4. — Radial expansion measure from the fits to the quadrilateral lobes plus torus model 
for the broken linear fit with an early fast shock velocity (va) transitioning to the later slower 
shock velocity (vb). The radii for the 3 image models are listed in Table [21 and the resulting 
radial expansion fits are given in Table [31 



